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SUMMARY 

An  experimental  investigation  of  the  laminar,  transitional,  and 
turbulent  heat  transfer  rates  over  a  conical  cylindrical  flared  body  is 
presented.  Regions  of  favorable,  zero,  and  adverse  pressure  gradient 
on  the  body  are  investigated.  The  experimental  results  are  compared 
with  the  theories  available  in  the  literature. 

The  modeL  chosen  for  this  investigation  is  a  cone -cylinder -flare 
configuration  consisting  of  a  20°  semi-vertex  conical  nose  portion 
smoothly  blended  by  a  shoulder  radius  into  a  long  cylindrical  body  and 
terminated  by  a  smooth  large  radius  flare. 

The  model  was  tested  at  a  free  stream  Mach  number  of  8  and  a 
Reynolds  number  of  1.6xlOB  to  0,3x1 06  per  inch  based  on  free  stream 
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conditions.  Various  stagnation-to-wall  temperature  ratios  were  obtained 
by  cooling  the  model  prior  to  the  test  with  Liquid  nitrogen.  The 
stagnation-to-waLL  temperature  ratios  ranged  from  10  to  3.3. 

The  theoretical  predictions  gave  good  resuLts  for  the  heat  transfer 
rates  in  the  Laminar  region,  and  fair  prediction  in  the  transitionaL  and 
turbulent  regimes  extending  over  the  shoulder  and  forward  portion  of  the 
cylindrical  body.  Over  the  aft  portion  of  the  cylinder  and  over  the  fLare 
the  predictions  are  only  qualitatively  correct,  and  underestimate  the 
heating  rate  by  a  factor  as  high  as  2.  ConverseLy,  the  "fLat  pLate  refer¬ 
ence  enthalpy"  method  is  found  more  closely  to  predict  the  heat  rates 
over  the  aft  portion  of  the  body,  but  increasingly  to  overestimate  the 
heating  rates  over  the  forward  portion  of  the  cy Linder. 

A  modified  equation  for  the  heat  transfer  coefficient  in  the 
transitionaL  and  fully  turbulent  region  based  on  the  F.P.R.E.  method 
is  then  presented.  This  method  gives  good  agreement  with  the  experi¬ 
mental  resuLts  presented  here. 

From  the  results  the  following  is  concluded:  Cooling  the  wall 
delayed  transition.  By  expanding  the  flow  rapidly  between  the  cone  and 
the  cylinder,  the  transition  ReynoLds  number  is  reached  very  rapidLy. 

By  making  a  smooth  transition  between  the  cylinder  and  the  fLare,  no 
separation  occurred  at  the  cylindrical  flare  junction.  The  transitional 
and  turbulent  heat  transfer  in  the  presence  of  an  adverse  pressure 
gradient  may  be  predicted  with  sufficient  accuracy  by  the  F.P.R.E. 
method. 
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SECTION  I 


INTRODUCTION 

In  recent  years*  with  the  progress  of  hypersonic  fLight,  the 
ability  reliably  to  predict  aerodynamic  heating  rates  has  become 
increasingly  important;  and  concomitant  with  the  developments  in 
rockets  and  ballistic  missiles,  considerable  attention  has  been  given  to 
the  study  of  the  boundary  layer  and  heat  transfer  rates  over  axially 
symmetric  bodies  in  high  speed  flow. 

Specifically,  the  aim  of  the  present  work  is  to  evaluate  the 
accuracy  with  which  several  of  the  more  recent  analytical  methods 
predict  boundary  layer  development  and  heat  transfer  rates  in  the  tran¬ 
sitional  and  turbulent  region  at  high  Mach  numbers.  Several  of  these 
evaluations  have  been  made  at  a  relatively  low  Mach  number,  and  it  is 
the  purposes  of  this  report  to  extend  these  evaluations  to  a  high  Mach 
number. 

For  the  laminar  flow  region  a  method  of  similar  flow  based  on 
the  numerical  solutions  of  references  1  and  2  is  employed  together 
with  the  similarity  transformations  of  reference  3.  For  the  transi¬ 
tional  and  fully  turbulent  regimes  a  modification  of  the  momentum 
integral  method  of  reference  4>  as  suggested  in  reference  5,  is  used 
together  with  the  Reynolds  anaLogy  corrected  empirically  for  nonunity 
Prandtl  number.  The  flat  plate  reference  enthalpy  method  has  also 
been  used  by  various  investigators  and  has  yielded  good  results  in 
some  cases. 
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After  formulation  of  the  methods,  solutions  are  obtained  for  a 
particular  cone-cylinder-flare  geometry,  and  finally  the  predicted 
results  are  compared  with  experiments.  However,  although  applied  to 
a  specific  body,  it  is  also  felt  that  the  results  are  of  general  interest 
in  ascertaining  the  applicability  of  the  methods.  For  example,  in  the 
present  work  a  region  of  particular  interest  is  the  junction  of  the 
conical  surface  with  the  cylindrical  body,  where,  because  of  the  radius 
of  curvature,  the  pressure  gradient  becomes  large,  and  thereby  pro¬ 
vides  a  critical  test  of  the  laminar  and  transitional  theory.  In  addition, 
the  compressing  flow  over  the  flare  portion  provides  for  examination  of 
heat  transfer  theories  under  the  influence  of  adverse  pressure  gradients. 

The  effect  of  a  cold  wall  is  studied,  and  its  effect  on  the  transi¬ 
tional  Reynolds  number  is  also  given. 

Some  experimental  heat  transfer  results  have  been  reported  in 
reference  6,  7,  and  8  from  free  flight  tests.  In  this  report  the  experi¬ 
mental  results  have  been  obtained  in  the  PIBAL  Mach  8  wind  tunnel. 

The  stagnation  pressure  ranged  from  600  to  100  psia,  and  the  stagnation 
temperature  was  approximately  2000°R.  The  Reynolds  number  based 
on  free  stream  conditions  ranged  from  1.6xl0*/into  to  .25xl0*/inch. 

This  work  has  been  initiated  by  Dr.  Antonio  Ferri,  and  the 
authors  wish  to  thank  him  for  his  suggestions  and  discussions  in  the  work 
reported  here.  Grateful  acknowledgment  is  due  to  Drs.  Paul  A.  Libby 
and  Robert  J.  Cresci  for  their  stimulating  discussion. 
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SECTION  II 


MODEL  DESIGN  AND  TESTING  PROCEDURE 

The  experimental  program  was  conducted  in  the  hypersonic 
facility  of  the  Polytechnic  Institute  of  Brooklyn  Aerodynamics  Laboratory. 
The  tests  were  conducted  at  a  Mach  number  of  8,  in  a  range  of  stagnation 
pressures  ranging  between  100  and  600  psia,  and  at  a  stagnation  temper¬ 
ature  of  approximately  2000°R.  The  stagnation  to  wall  temperature  of  the 
model  was  varied  from  Tq/Tw  from  3.2  to  10  by  cooling  the  skin  temper¬ 
ature  of  the  model  with  liquid  nitrogen  prior  to  the  test.  The  range  of 
free  stream  Reynolds  number  per  inch  for  the  tests  presented  here  is 
1.6xl06/in.  x  0.25x1 0E/in. 

A  schematic  drawing  of  the  model  tested  is  shown  in  Fig.  1.  The 
model  consists  of  a  pointed  conical  portion  having  a  half  angle  of  20 
degrees.  This  conical  portion  terminates  at  an  S  of  4.7  (where  S  is  S/Ro 
and  Rq  is  a  reference  Length  taken  to  be  1").  The  conical  portion  is 
blended  within  a  cylindrical  portion  by  means  of  a  segment  of  a  circuLar 
arc.  The  equation  representing  the  geometry  of  this  circular  arc  in 
terms  of  r  is  given  as 

4.696  <  S  <  6.616 

7  =  -  3.563  +  5.5  cos  [68.92-10. 42S] 

The  circular  arc  terminates  with  the  cylinder  at  S=6.616.  The 
cylindrical  portion  continues  then  to  an  S  of  23.99.  At  this  point  the 
cylindrical  body  flares  out  from  a  diameter  of  3.875"  to  a  diameter  of 
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5.95"  in  a  distance  of  12".  Similarly  in  this  case  the  flare  is  also 
expressed  by  a  segment  of  a  circle  and  is  given  in  terms  of  r  by  the 
following  equation: 

23.99  <S<  36.19 

r  =  71.854  -  69.916  cos[19.66  -  0.82S  ]  . 

The  model  was  constructed  of  304  stainless  steel  having  a  wall 
thickness  of  0.050".  Table  II  presents  the  location  of  the  thermo¬ 
couples  and  the  pressure  taps.  The  pressure  taps  and  thermocouples 
were  located  at  two  planes  at  a  radial  angle  of  30°.  A  copper  tube, 
having  a  one-half  inch  diameter  and  having  several  holes  1/32"  in 
diameter,  was  placed  at  the  center  line  of  the  model  along  all  its 
length  for  the  purpose  of  cooling  the  model  with  Liquid  nitrogen  prior 
to  the  test. 

The  heat  transfer  to  the  modeL  was  determined  by  the  transient 
technique.  The  model  is  initially  at  constant  temperature;  then  the 
initial  heat  pulse  was  obtained  by  starting  the  flow  suddenly.  The 
derivative  of  the  surface  temperature  with  respect  to  time  was 
measured  at  0.2  seconds  from  starting  and  used  in  determining  the 
heat  transfer  to  the  surface  of  the  model.  For  this  condition  the  temper¬ 
ature  of  the  model  was  almost  uniform,  therefore  lateral  conduction 
could  be  neglected.  In  the  case  where  a  temperature  less  than  the 
ambient  temperature  was  desired  prior  to  the  test,  the  following  pro¬ 
cedure  was  followed.  The  tunnel  was  evacuated  to  a  vacuum  of  29 
inches  of  mercury,  and  the  copper  tube  within  the  model  was  connected 
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to  a  liquid  nitrogen  bottle.  Cooling  of  the  model  continued  until  a 
uniform  temperature  on  the  model  was  reached.  By  maintaining  a 
vacuum  in  the  wind  tunnel,  no  frosting  occurred  within  or  at  the 
surfs  ce  of  the  model. 

The  data  are  presented  in  terms  of  a  Nusselt  number,  defined 
by 

NNU  *  V<Cp'.eRoA..<h..-hw' 
and  by  a  Reynolds  number  parameter 

N_,  =  N_  </p"  /P  h  =  N„cp  * 

R  R  s'  se  se  R^se 

N_  =  p  Jh  R  /n 
R  Kse  se  o'  se 

cp  =  P  /p  h 
^se  se' Kse  se 

Four  test  conditions  are  presented  in  the  report, 


Run  No, 

P  (psia) 

00 

T 

w 

1 

590 

520 

2 

200 

520 

3 

100 

520 

4 

600 

180 

All  the  tests  were  performed  at  a  stagnation  temperature  of 
approximately  1800°R. 
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SECTION  III 


PRESSURE  DISTRIBUTION  AND  HEAT  TRANSFER  RESULTS 


The  pressure  distribution  resuLts  are  presented  in  Fig.  2  in  the 
form  of  P  /P  versus  "§,  where  P  is  the  stagnation  pressure  behind  the 
conical  shock.  The  Mach  number  on  the  conicaL  portion  is  Me=4.24, 
which  agrees  well  with  the  Kopal  value.  A  Mach  number  of  6.80  is 
reached  after  the  expansion;  this  value  is  somewhat  higher  than  the  value 
obtained  with  a  standard  Prandtl-Meyer  expansion.  This  is  due  to  the 
overexpansion  which  occurs  at  the  beginning  of  the  cylinder.  At  the  end 
of  the  flare  the  flow  is  compressed  from  a  Mach  number  of  6.8  to  a 
Mach  number  of  5.84.  Included  in  the  figure  is  the  comparison  of  the 
pressure  distributions  on  the  flare  with  the  tangent  cone  method  (Table  I). 

The  experimental  heat  transfer  data  are  presented  in  three 
different  forms  in  order  to  distinguish  the  effects  of  laminar,  transi¬ 
tional,  and  turbulent  heat  transfer.  The  results  are  presented  in  terms 
of  the  Nusselt  and  Reynolds  numbers  due  to  their  insensitivity  with 
respect  to  enthalpy  ratios. 

_  \ 

It  is  well-known  that  for  laminar  flow  N..  varies  with  N„2 

Nu  R 

~  4/5 

while  for  turbulent  flow  the  variation  is  N.T  ~  N_  '  . 

Nu  R 

In  order  to  indicate  the  laminar  part  of  the  variation  the  heat 
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transfer  data  is  presented  in  Figs.  3-5  in  the  forms  of  N^ju/N^2 
versus  S.  Where 


(h  -h  )K ' 
'  aw  w 
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(2) 


and 


py 


S  =  -5-  (R  =1) 
R_  ° 


It  is  seen  that  by  presenting  the  data  in  this  form  both  the 
Laminar  part  and  transitional  point  on  the  body  may  easily  be  indicated. 
Figs.  6-9  present  the  variation  of  N^u  with  for  laminar  transitional 
and  turbulent  heat  transfer  results  for  the  full  body  length.  For  com¬ 
parison  with  the  various  theories  the  experimental  results  are  again 
presented  in  Figs.  10  and  11  for  test  no.  1  in  the  form  of  Nj^u  and 
N'  versus  S,  where 


N 


Nu 


q  C  R 
w  p  o 
rse 

k  (h  -h  ) 
se  se  w' 


(3) 


SECTION  IV 

THEORETICAL  ANALYSIS 

The  more  commonly  employed  method  for  computing  the  Laminar 
boundary  layer  is  to  use  the  transformation  of  reference  3  for  similar 
solutions,  with  the  numerical  solution  of  reference  1.  The  above  method 
yields  good  results  when  the  pressure  gradient  as  defined  in  the 
equation  below,  does  not  vary  significantly  along  the  body: 
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(4) 


where 


r. 


dM  h 

-  <Tr'  ^.e. 
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dS 


■5 

pLees 


2S 
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du 
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dS 


(5) 


For  Larger  vaiues  of  "j?(T?>2)  the  numerical  solutions  of  refer¬ 
ence  2  have  to  be  used.  The  initial  estimates  of  Laminar  heat  transfer 
were  based  on  the  assumption  that  "B=  0. 

In  terms  of  the  nondimen sionaL  quantities  given  in  reference  5, 
the  theoretical  results  for  the  Laminar  heat  transfer  of  reference  3  give 


N-. 

Nu 

“  I 
N  2 
NR 


J2~  -j  °w 


g 


ppur 


3 

,3 


2  '-Sw  [JSpMar'dS]*  VJ 
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for  the  case  "5=0,  the  similar  solutions  of  reference  l  give 


gw 

(_rL_)  =  0.470 

1-g 

°w 


(7) 


^Nu  _  0.67  p  |i  u  r _  03  (8) 

V  2  cp  * 

o 

Over  the  conical  nose  portion  of  the  model  this  result  reduced  to 
the  simple  form 
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(9) 


Mi  /r  (i^)1 
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From  these  computations  the  predicted  Laminar  heating  ratio 
for  the  present  body  are  presented  in  Fig.  10.  The  agreement  between 
these  results  and  the  experiments  are  seen  to  be  good. 

The  laminar  heat  transfer  over  the  conical  portion  of  the  modeL 
may  aLso  be  calculated  in  terms  of  the  reference  enthalpy  .method.  In 
terms  of  the  notation  presented  in  reference  5,  the  Laminar  heat 
transfer  is 
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Nu 
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Nu 
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=  0.576  a'3  (N 
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(flat  pLate) 

(cone) 
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(prime  quantities  indicate  properties  evaluated  at  the  reference 
conditions). 

The  resuLts  of  these  analyses  are  presented  in  Figs.  2-9  for 
the  tests  reported  here. 

A.  Laminar  momentum  thickness  and  the  transition  point 

At  the  outset  of  this  work,  it  was  known  from  experimental 

observations  that  transition  from  laminar  to  turbulent  flow  occurred 

approximately  at  a  Local  Reynolds  number,  based  on  momentum 

a 

thickness,  equal  to  approximately  600.  To  establish  the  location  of 


9 


the  transition  point  for  the  theoretical  anaLysis  it  became  necessary  then 
to  calculate  the  development  of  the  laminar  momentum  thickness  over 
the  body. 


l.  Solution  by  exact  similarity 

From  the  definition  of  the  momentum  defect  thickness 


-J 


oo 


e*-  (l-JL)dy 


n  p  u  u 

°  e  e  e 


(12) 


and  the  locaL  Reynolds  number 


P  u 

"Re  -  tP  • 


(13) 


the  transformation  of  Lees  results  in 


e  .  138.  *-  e 


p  u  r 
e  e  o 


(14) 


and 


N  fl  =  l&i  e 

u  r 
e  o 


(15) 


where  ©  is  the  transformed  momentum  thickness 


0  = 


f '(1-f  #)dTi 


(16) 


Employing  the  method  of  similar  flows,  the  value  of  9  is  com¬ 
puted  from  the  solution  of  the  velocity  profile  f / ,  and  hence  0  is 
functionally  dependent  on  the  similarity  parameters  *6  and  gw> 
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In  the  present  work,  for  the  500°R  wall  and  1800°R  stagnation 
temperatures  upon  which  the  analytical  treatment  is  based,  the 
enthalpy  ratio  gives 


g  =  0.2772, 

°w  ' 

and  as  computed  from  the  similar  solutions  of  references  1  and  2, 

Fig.  12  presents  the  corresponding  values  of  0  as  a  function  of  the 
pressure  gradient  similarity  parameter  "5. 

Over  the  conical  portion  of  the  body,  the  relation  between  body 
radius  and  the  surface  coordinate  is  given  explicitly  by 

r  =  S  sin  6 
o  c 

Hence,  the  transformed  surface  coordinate  S  becomes 

S'K.V1.  T  ““’‘c  (17) 

and  by  Eqs.  (14)  and  (15),  the  physical  momentum  thickness  and  local 
Reynolds  number  are  given  by 


M, 


P  u 
e  e 


S>*  ® 


and 


p  u  1 

sr  e 


( I  8'. 


(19) 


where  the  transformed  momentum  thickness  0  is  evaluated  from  Fig.  3 
for  "5=0. 
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Over  the  shoulder  region,  however,  the  properties  of  the  external 


flow  are  dependent  on  the  assumed  Prandtl-Meyer  expansion.  Further, 
a  simple  expression  for  rQ  as  a  function  of  S  is  not  easily  written.  Hence, 
in  this  region  S,  and  thus  0  and  N^g  ,  had  to  be  computed  by  numerical 
integration. 

For  this  calculation,  in  keeping  with  the  assumption  of  exact 
similarity,  the  value  of  ?  was  supposed  constant  all  along  the  body  surface. 
Hence,  while  in  fact  the  value  of  ?  varied  locally,  the  transformed 
momentum  thickness  ®  was  evaluted  (as  for  the  cone)  at  ?=  0. 

The  caLcuLated  distribution  of  N^g  over  the  forward  portion  of 
the  body,  based  on  the  presently  prescribed  external  flow,  is  shown  in 
Fig.  13a. 

According  to  the  transition  criteria,  600,  Fig.  13a  shows 

the  predicted  transition  point  to  occur  on  the  shoulder  at  S=6.07'!.  This 
location  on  the  body  is  illustrated  in  Fig.  1. 


2.  Solution  by  piecewise  application  of  exact  similarity 

As  mentioned  in  the  above  analysis,  the  flow  expansion  over 
the  cone-cylinder  shouLder  region  results  in  a  Large  variation  in  the 
Local  value  of  the  pressure  gradient  parameter  ?.  Since  the  previous 
analysis  predicted  the  transition  point  to  Lie  in  this  region,  it  became 
of  interest  to  investigate  the  effect  of  f  more  rigorously. 

The  approach  often  employed  to  account  for  the  effect  of 
variations  in?  is  the  method  of  simple  iocaL  similarity.  However, 
in  that  method,  since  ?  depends  only  on  the  inviscid  flow  field,  there  is 
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no  accounting  for  the  previous  history  of  boundary  Layer  development. 

To  account  in  part  for  previous  boundary  Layer  history,  it  was 
decided  to  employ  the  method  of  piecewise  exact  similarity.  This 
method  insures  continuity  in  the  distribution  of  boundary  Layer  momentum 
thickness,  and  is  based  upon  an  approach  first  suggested  for  incom¬ 
pressible  fLows  by  A.M.O.  Smith. ^  The  present  development  follows 
cLoseLy  that  of  reference  2,  except  that  in  that  case  heat  transfer  was 
of  primary  concern,  and  the  development  was  predicted  on  continuity 
in  the  distribution  of  the  convection  thickness. 

The  basic  assumption  underlying  the  present  piecewise 
application  of  exact  similarity  is  as  folLowsr 

"An  exact  similar  flow,  not  necessarily  originating 
at  the  same  point  as  the  actual  fLow,  may  be  appLied 
within  a  finite  intervaL  of  the  body  surface,  and  the 
effect  of  previous  boundary  Layer  history  may  be 
introduced  by  requiring  continuity  in  the  physical 
momentum  thickness," 

Specifically,  let  us  consider  a  surface  intervaL  [1]~[2].  Within  this 
intervaL  the  external  veLocity  distribution  is  known  from  the  potential 
flow  field.  Now,  the  physicaL  momentum  thickness  and  the  external 
veLocity  at  the  start  of  the  interval  may  be  matched  by  a  simiLar  flow 
of  arbitrary  origin,  and  the  arbitrariness  of  the  origin  may  then  be 
removed  by  requiring  that  the  simiLar  fLow  also  satisfy  the  external 
velocity  at  the  end  of  the  intervaL.  This  simiLar  flow  may  then  be  used 
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to  predict  the  physical,  momentum  thickness  at  the  end  of  the  interval. 

Starting  with  a  modified  form  of  the  transformation  of  Lees,  for 
axisymmetric  bodies 


P  u 
Ke  ' 


u 

*  *  J  ro  T"  dy 

v  2(S-S_)  O  Pe 


and 


S  =  [  pelieu 


r  fldS  , 

e  o  ’ 


(20) 


(21) 


where  SQ  will  be  the  effective  origin  of  the  similar  flow. 
The  expression  for  IS"  becomes  in  this  case 

2(S-S  )  dM 


M 


z  > 
dS 


(22) 


and  hence,  what  has  been  said  previously  in  connection  with  exact 
similar  solutions,  still  applies. 

Now,  consistent  with  the  assumption  of  exact  similarity 
within  a  finite  interval,  ~B  and  SQ  are  required  to  be  constant.  Hence, 
upon  integrating  Eq.  (22)  within  an  interval  [l]-[2] 


S  -S  1 

[1 

Si"So 


(23) 


ALso,  employing  the  modified  similarity  transformation,  the  trans¬ 
formed  momentum  thickness  ©  must  also  be  constant  within  the  interval, 
and  thus  for  the  end  points  of  the  interval 
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(24) 


and 
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Finally,  upon  manipulation  of  these  relations,  the  working  equations  for 
the  interval  are  obtained  as 


M 


M 
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•  si 


■  {i  +  afr-S.)  [?■■-]} 
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and 


<peVo>i  “v* 
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<PelVo>  Me 

3  1 


>  Bi 


(27) 


The  first  of  these  equations  is  soLved  by  a  trial-and-error  procedure 
for  H,  making  use  of  the  dependence  of  0  on  ;  and  after  obtaining  TS, 
the  physical  momentum  thickness,  0a ,  at  the  end  of  the  interval  is 
found  from  the  second  equation. 

For  the  present  work,  in  the  cone-cylinder  shoulder  region, 
application  of  this  method  of  "Piecewise  Exact  Similarity"  resulted  in 
the  schedule  of  similarity  parameter  7  shown  in  Fig.  14.  The  corre¬ 
sponding  distributions  of  the  Reynolds  number  N^g  has  been  super¬ 
imposed  on  the  plot  of  Fig.  13a. 

From  Fig.  13a,  again  employing  the  transition  criteria 
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N^g  -  600,  the  location  of  the  transition  point  as  predicted  by  the  method 
of  Piecewise  Exact  Similarity  is  found  to  occur  at  the  station  S=5.60 
inches;  for  reference,  this  location  has  been  indicated  on  the  configu¬ 
ration  drawing  of  Fig.  1. 

Comparing  the  reLative  locations  of  the  transition  point  as 
predicted  by  the  method  of  exact  similarity  and  the  piecewise  method 
of  exact  similarity  respectively,  the  difference  is  seen  to  be  smaLl. 

The  sLight  difference  in  the  predicted  Location  of  the  transition 
point  was  considered  to  be  welL  within  the  accuracy  of  the  N^g  —  600 
transition  criteria  employed.  Hence,  it  was  concluded  that  the  rigor 
of  the  method  of  piecewise  solution  was  unnecessary,  and  from  this 
stage  on,  the  method  of  exact  similarity  was  considered  sufficiently 
accurate  for  all  calculations  associated  with  the  Laminar  boundary 
Layer. 

The  distribution  of  Nj^g  in  Fig.  13b  has  been  calculated  for 
test  no.  1.  In  order  to  obtain  the  distribution  of  N^g  over  the  body 
for  the  other  test  conditions  the  following  equation  is  used: 


Nr0s  8°°a  i 

=  <— -) 


p 

i  soo , 


(28) 


By  applying  the  above  equation  at  each  point  of  the  body,  the 
Laminar  distribution  of  N^g  corresponding  to  any  other  free  stream 
stagnation  pressure  p  may  be  found. 

8Cj| 


16 


The  effect  of  free  stream  stagnation  pressure  on  the  transitional 


point  for  N^g  =600  has  been  calculated  from  the  above  equation  and  is 
presented  in  Fig.  15. 

From  Figs.  13b  and  15  it  is  seen  that  the  transitional  point 

persists  for  a  larger  range  of  stagnation  pressures  over  the  expansion 

section  of  the  body  than  on  the  cylindrical  section  of  the  body. 

It  was  of  considerable  interest  to  note  that  by  expanding  the  flow 

rapidly  around  the  cone  N_  a  increases  very  rapidly,  and  reached  a  value 

x\  y 

of  600  at  S  equal  to  5.6.  On  the  other  hand,  if  the  cone  extended  further, 
the  transitional  point  would  occur  at  an  S  of  13,4  (this  is  consistent  with 
the  results  presented  in  reference  9).  These  results  are  shown  in  Fig.  13b. 

In  the  case  of  run  4,  because  of  the  cold  wall  conditions,  transition 
occurs  at  an  N^g  of  1100,  Cooling  the  wall  thus  has  an  effect  on  delaying 
transition,  as  has  been  observed  by  previous  investigation.^ 

The  following  table  gives  a  summary  of  the  transitionaL  point  on 
the  bodies  for  the  tests  reported  here: 


Test  No. 

S 

P 

o 

DO 

T  /T 
os'  w 

NR0 

1 

5.6 

590 

3.3 

600 

2 

15 

200 

3.3 

680 

3 

23 

100 

3.3 

450 

4 

24 

590 

10 

1100 

B.  The  transitional  and  turbulent  wall  friction  coefficient 

To  evaluate  the  wall  friction  coefficient  in  the  transitionaL  and 
fully  turbulent  flow  regimes,  the  friction  law  modified  in  reference  5  was 
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adopted,  i.e.r 


0.013 


B 


N. 


nl 


R0 


(29) 


4 

This  Law  is  simiLar  in  form  to  that  developed  by  Persh  through 
anaLogy  with  incompressible  flow,  but  has  been  modified  to  include  the 
law  suggested  by  Vaglio-Laurin^  for  the  fully  developed  turbulent 
fLow  of  a  compressibLe  medium. 

The  constant  B  of  this  law  is  evaluated  from  the  laminar  flow. 


such  that  the  friction  coefficient  —  remains  continuous  at  the  laminar 

2 

to  turbulent  transition  point,  i.e,: 


where 


and 


and  therefore 
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In  the  analysis  of  reference  5  the  exponent  n  has  been  chosen  to  be  2, 
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Recently,  reference  13  has  indicated  that  the  use  of  n=l  leads  to  better 
agreement  for  the  transitional  heat  transfer  in  a  favorable  pressure 
gradient  region.  Therefore  in  the  following  analysis,  in  the  favorable 
pressure  gradient  region,  theoretical  calculations  were  carried  out  for 
n=l  and  2  for  one  of  the  test  conditions,  and  are  presented  in  Table  II. 

It  is  seen  that,  in  this  case,  the  variation  of  n  on  the  transitional  heat 
transfer  has  LittLe  effect.  The  reason  for  this  is  that  the  pressure 
gradient  (?)  effect  in  this  case  is  much  more  predominant  than  the 
effect  of  n.  In  the  case  of  reference  13  the  Mach  number  and  pressure 
gradient  were  low  in  comparison  to  the  results  presented  here. 

C.  The  transitionaL  and  turbulent  boundary  layer 

1,  Modification  of  the  method  of  Persh 

For  prediction  of  the  transitional  and  turbulent  boundary  layer 

behavior,  with  inclusion  of  the  effect  of  compressibility,  body  geometry, 

and  pressure  gradient,  the  present  work  employed  a  calculation  scheme 

4 

based  on  a  modification  of  the  method  introduced  by  Persh. 

According  to  the  method  of  Persh,  a  stepwise  solution  of  the 
boundary  layer  momentum  equation 


do 

dS 


C,  •*-+2du  .  dp  ,  dr 

_J  .  0[° - «  ♦  i - i  +  i - °] 

2  u  dS  p  dS  r  dS 


(34) 


is  obtained  through  use  of  boundary  layer  form  factors,  which  are  based 
on  a  correlation  of  incompressible  data. 

In  the  present  work  however,  following  the  suggestion  of 
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reference  5,  the  form  factor  was  taken  as 


With  this  assumption  thep,  the  momentum  equation  reduced  to  the 
modified  form 


C,  ,  du  ,  dp^  ,  dr 

M  =  -i  _  0[i _ -  +  - _ 11  +  J _ ^1 

dS  2  u  dS  p  dS  r  dS 
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(35) 


and  combining  this  expression  with  the  wall  friction  law,  the  Local 
wall  friction  coefficient  and  the  boundary  layer  development  may  be 
obtained  by  stepwise  numerical  solution. 

An  approximate  method  for  soLution  is  discussed  in  detail 
below.  However,  as  stated  in  the  general  remarks,  boundary  Layer 
solutions  in  the  present  work  were  obtained,  together  with  heat 
transfer  rate  predictions,  on  the  Bendix  automatic  computer.  From 
these  computations,  the  distribution  of  Reynolds  number,  , 

corresponding  to  each  of  the  pressure  fields,  is  shown  in  Fig.  13b:  and 
in  addition,  the  tabulated  vaLues  are  presented  together  with  the  heat 
transfer  resuLts  in  TabLe  II. 


*ALthough  this  was  the  case  assumed  for  bLunt  bodies  it  was  of 
interest  to  see  the  resuLts  of  such  an  analysis  when  applied  to 
this  type  of  body. 


20 


D.  Transitional  and  turbulent  heat  transfer 
1.  The  Reynolds  analogy 

To  compute  heat  transfer  rates  in  the  regions  of  transitional  and 
fully  turbulent  flow,  the  Reynolds  analogy,  corrected  empirically  for 
the  effect  of  nonunity  Prandtl  number,  may  be  written 


g=CT*(h  -h  ) p  u  -~ 
°w  '  aw  wee  2 


(36) 


where  the  adiabatic  wall  enthalpy  is  defined  in  terms  of  the  recovery 
factor 


=  a3 


and  the  wall  friction  coefficient  is  evaluated  from  the  boundary  layer 
analysis. 

Framing  this  result  in  terms  of  the  Nusselt  and  ReynoLds 
numbers  previously  defined,  we  write 
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and  finally 
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E.  Heat  transfer  expressed  in  reference  enthalpy  form 

For  real  gases,  although  properties  across  the  boundary  layer 
are  known  to  vary  with  temperature,  the  simple  assumption  of  a 
constant  property  fluid  has  been  found  to  describe  actual  heat  transfer 
processes  with  good  accuracy,  providing  the  properties  are  introduced 
at  a  suitable  reference-enthalpy  state.  For  this  reason  then,  in 
reducing  the  experimental  data  to  nondimensional  form,  local  Nusselt 
and  Reynolds  numbers  were  defined  as 


and 


N '  = 

Nu  k'(h 


g  C'S 

°W  p 


-h  ) 
aw  w 


(39) 


(40) 


where  (for  example,  see  reference  5),  the  primed  quantities  are 
evaluated  at  the  state  corresponding  to  the  reference  enthalpy: 

I  i. 

h'  =  0.5h  +  0.22a3  h  f  h  (0.50-0. 22a3)  (41) 

w  se  e' 


1.  Conversion  of  the  analytical  results 

In  order  to  compare  the  analytically  predicted  heat  transfer 
results  to  the  results  of  the  experiment,  it  was  necessary  to  convert 
the  predicted  results  to  the  reference-enthalpy  form  defined  above. 
For  this  purpose,  the  Nusselt  number,  N^u  was  rewritten  as 
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(42) 


The  Local  values  of  the  quantity  (M  thus  provided  for  conversion 

from  the  Nusselt  number,  N^u  used  in  the  computer  program,  to  the 
reference-enthalpy  Nusselt  number,  N^u  in  which  the  experimental 
results  were  expressed. 

2.  Comparison  of  predicted  and  experimental  heating  rates 

From  the  analytical  results  obtained,  the  predictions  of  laminar, 
transitional,  and  turbulent  heating  rates  over  the  cone-cylinder-flare 
configuration  are  presented  in  Fig.  11,  in  terms  of  the  reference- 
enthalpy  Nusselt  number  N^u.  The  predictions  are  presented  for  the 
experimental  pressure  field  computations,  and  from  Table  IV  the 
results  of  the  experimental  investigation  are  shown  for  comparison. 

From  Fig.  II ,  it  is  seen  that  in  the  laminar  flow  regions,  over 
the  conical  nose  and  the  beginning  of  the  cone -cylinder  shoulder,  the 

t 

analytical  predictions  match  the  experimental  values  quite  well.  Also, 
in  the  transitional  and  turbulent  regimes  on  the  shoulder  and  over  the 
forward  part  of  the  cylinder,  the  correlation  is  fairly  good.  However, 
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over  the  after  portion  of  the  cylinder  and  over  the  flare,  the  predicted 
heating  rates  are  only  qualitatively  correct,  and  increasingly  under¬ 
estimate  the  actual  turbulent  heating  rate, 

3.  The  F.P.R.E.  method 

In  the  investigation  of  reference  5,  based  on  the  heat  transfer 
over  a  spherically  capped  cone,  the  flat  plate  reference  enthalpy 
•  (F.P.R.E.)  method  was  found  to  give  quite  accurate  predictions  of 
fully  developed  heating  rates.  It  was  thus  of  interest  to  determine 
whether  this  simple  method  would  result  in  a  satisfactory  prediction 
of  heating  rates  in  the  turbulent  regime  of  the  present  cone-cylinder- 
flare  configuration. 

On  the  basis  of  the  local  reference  enthalpy  Nusselt  and 
Reynolds  numbers  used  herein,  the  basic  formula  for  the  F.P.R.E. 
method  is: 

X  4 

N^u  =  0.030  c^N^)*  (43) 

The  results  of  the  computation  for  the  present  problem  have  been 
added  to  the  plot  of  Fig.  11,  and  although  the  heating  rate  is  more 
closely  predicted  over  the  aft  portion  of  the  body,  the  F.P.R.E. 
method  is  seen  consistently  to  overestimate  the  heating  rate  over  the 
forward  portion  of  the  cylinder. 

4.  T  ransition  point  adjustment 

From  inspection  of  the  experimental  heating  rate  distribution 
(see  Fig.  3a),  a  sharp  increase  is  seen  to  occur  at  approximately  the 
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station  S=8.4  inches.  Since  this  behavior  is  generally  characteristic  of 
transition  to  turbulent  flow,  it  was  suspected  that  the  previous  conclusions 
concerning  the  location  of  the  laminar  to  turbuLent  transition  point  might 
have  been  in  error.  Accordingly,  to  investigate  this  possibility,  the 
transition  point  was  empirically  adjusted  to  this  location,  and  the  com¬ 
plete  turbuLent  boundary  Layer  and  heat  transfer  calculation  was  repeated, 
using  the  BendixG-15  computer  routine.  The  results  of  this  computation 
are  added  to  the  plot  of  Fig.  11,  and  as  seen,  no  significant  change  in 
the  heat  transfer  rate  prediction  is  obtained.  In  fact,  the  resulting  pre¬ 
dictions  are  slightly  lower.  The  tabulated  results  are  presented  in 
Table  HI. 

F.  Present  analysis 

It  has  been  shown  in  the  previous  section  that  the  modification  of 
the  method  of  Persh,  as  presented  by  reference  5,  yielded  an  unsatis¬ 
factory  result.  This  may  be  attributed  to  the  fact  that  the  method  of 
reference  5  applies  only  to  blunt  bodies  which  have  a  Low  surface  Mach 
number,  as  in  the  case  of  the  Laminar  theory  of  Lees.  In  this  theory  a 

form  factor  of  H  =  -  1  was  assumed,  consistent  with  the  cold  wall 
c 

assumption  of  Lees.  One  can  envisage  a  negative  displacement  thickness 
on  a  blunt  body;  however,  for  the  body  under  investigation,  this  is  not 
possible,  since  the  Local  surface  Mach  number  is  fairly  high.  Based  on 
the  above  conclusions,  and  the  fact  that  the  F.P.R.E.  method  yield  good 

results  in  the  fully  developed  turbulent  region  on  the  fLare,  a  method 

■  -■  -■■■■■■■■■-■■  -  ■  -  ■  ■  ■  1  —  ■  —  — —  ~  "  ~  '  ’  ■  1 

Extension  of  the  laminar  calculation  to  S=8.4"  gave  (NRSltrans” 802,  for 
start  of  this  calculation. 
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similar  to  method  5  was  adopted  for  the  transitional  region,  onLy  in  this 
case  the  F.P.R.E.  method  was  adopted  throughout  the  transitional  and 
turbulent  calculations.  This  will  ensure  that  the  value  of  the  fully 
developed  turbulent  boundary  layer  will  match  the  value  as  calculated  by 
the  F.P.R.E.  method.  Furthermore,  the  present  method  requires  no 
integration  of  the  momentum  equation  and  may  be  calculated  rapidly  to 
yield  the  transitional  and  turbulent  heat  transfer  on  the  body.  In  terms 
of  the  reference  enthalpy  method,  the  following  procedure  is  foLLowed. 

The  skin  friction  in  terms  of  the  reference  enthalpy  method  is 
given  as 


_ f  _  0.332  (Laminar  Flowj 
2  <NR>5 

Q 

—  =  (.0296)/n'  )°*2  (Turbulent  Flow) 

b  ^ 


(44) 


(45) 


To  evaLuate  the  wall  friction  coefficient  in  the  transitionaL  and 
fully  turbulent  flow  regimes,  the  skin  friction  low  of  reference  5  may 
be  modified  and  written  in  the  following  way: 


.0296 


B. 


K 


0.2 


N. 


/  a 


(46) 


where  n=l  for  pressure  gradient  and  n=2  for  zero  pressure  gradient. 
Eq.  (46)  may  be  written  in  terms  of 
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<47) 


NNu  =  °.°30o*(Nr);- 

1NR 

In  this  case  the  constant  B  is  evaluated  such  that  the  heat  transfer 
NNu  8iven  by  Eq.  (46)  is  continuous  at  the  laminar  to  transitional 

1  4 

B  =  N^[.03o3N^5-N^u]  (48) 

In  the  case  of  test  no,  1  transition  occurs  at  the  beginning  of  the 
cylinder. 

(Nj^)  Transition  =  4x10* 

(I*'  )  Transition  =  50 

Therefore  B  may  be  determined  for  Eq.  (48) 

B  =  350x10* 

The  result  of  such  an  analysis  is  included  in  Fig.  11.  It  is  seen 
that  very  good  agreement  is  obtained  for  this  method  when  compared  to 
the  experimental  results. 


SECTION  V 
CONCLUSION 

An  analytical  investigation  has  been  conducted  for  prediction 
of  the  boundary  layer  behavior  and  heating  rates  over  a  cone-cylinder 
flare  configuration  at  a  Mach  number  ahead  of  the  body  equal  to  8.0, 


27 


and  with  a  free  stream  to  model  wall  enthalpy  ratio  equal  of  3.6  and  10. 

Based  on  a  comparison  of  the  theoretical  heat  transfer  pre¬ 
dictions  with  experimental  measurements,  it  is  found  that  excellent 
results  are  obtained  in  the  laminar  region  by  using  the  transformation 
of  Lees  and  the  method  of  similar  solutions.  Even  in  the  shoulder 
region  joining  the  cone  and  cylinder,  the  predictions  are  quite  good  in 
spite  of  the  presence  of  severe  favorable  pressure  gradients. 

In  the  turbulent  flow  region  extending  over  the  cylinder  and 
flare,  however,  attempts  at  predicting  the  heat  transfer  rate,  based 
on  the  boundary  layer  behavior  obtained  from  a  modification  of  the 
momentum  integral  method  of  Persh,  were  rather  discouraging. 
Although  the  results  obtained  were  found  to  be  qualitatively  correct, 
they  consistently  underestimated  the  actual  heating  rates;  the  agree¬ 
ment  became  poorest  over  the  after  portion  of  the  cylinder  and  over 
the  flare.  Hence,  for  practical  prediction  of  the  transitional  and 
turbulent  heating  rates,  the  modified  momentum  integral  method  did 
not  yield  satisfactory  results. 

A  modified  equation  for  the  heat  transfer  coefficient  in  the 
transitional  and  fulLy  turbulent  region  based  on  the  F.P.R.  E.  method 
is  presented.  This  method  gave  good  predictions  of  the  heat  transfer 
over  the  entire  range  of  transition  and  fully  developed  turbulent 
region. 

From  these  results  it  is  concluded  that  cooling  the  wall 
delayed  transition.  By  expanding  the  flow  rapidly  between  the  cone 
and  the  cylinder,  the  transition  Reynolds  number  is  reached  very 
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rapidly.  The  transitional  and  turbulent  heat  transfer  in  the  presence  of 
an  adverse  pressure  gradient  may  be  predicted  with  sufficient  accuracy 
by  the  F.P.R.E.  method. 
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TABLE  I 


THE  ANALYTICAL  FLOW  FIELD 


M„  =  8.0  ,  Pc  -  600  psia 

00  5 CD 

Tc  =  1800°R,  T  =  500°R 
&oo  w 


A.  Conical  Nose: 

0  <  S  <  4.695,  P  =  0.4745x10 


B.  Cone-Cylinder  Shoulder: 


s 

PxlO3 

4.695 

0.4745 

4.887 

0.3815 

5.079 

0.3045 

5.271 

0,2406 

5.463 

0.1886 

5.655 

0.1461 

5.847 

0.1121 

6.039 

0.0848 

6.231 

0.0635 

6.423 

0.0468 

6.615 

0.0340 

C.  Cylinder: 


6.615  <S  <  23.99,  F  =  0,0340x10'* 


D.  Flare: 


S  PxlO3 


23.99  0.0340 

36.19  0.108 
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TABLE  II 


PREDICTION  OF  TRANSITIONAL  AND  TURBULENT 
BOUNDARY  LAYER  DEVELOPMENT  AND  HEAT  TRANSFER  RATES 

(Theoretical  Transition  Point,  S=6.07")  vtsX 


n=l 


n=2 


S 

NR0 

N.. 

Nil 

NR0 

N,t 
Nu  j 

6.12 

615.5 

651.9 

615.6 

6.560 

6.22 

629.7 

599.0 

629.8 

6.058 

6.32 

643.7 

545.7 

643.8 

5.545 

6.42 

655.6 

507.7 

655.8 

5.178 

6.52 

666.5 

477.4 

666.8 

4.884 

6.62 

677.2 

451.4 

677.5 

4.631 

7.12 

746.0 

370.8 

746.8 

3.858 

8.12 

822.4 

283.6 

824.3 

2.981 

9.12 

867.6 

245.1 

870.5 

2.587 

10.12 

899.3 

226.0 

903.2 

2.391 

11.12 

922.0 

217.8 

926.8 

2.307 

12.12 

940.0 

215.1 

945.7 

2.281 

13.12 

958.7 

210.4 

965.4 

2.233 

14.12 

974.9 

208.9 

982.5 

2.218 

15.12 

991,2 

206.8 

999.7 

2.198 

16.12 

100 

6 

206.4 

1015.4 

2.194 

17.12 

102 

0 

206.3 

1030.6 

2.194 

18.12 

103 

2 

208.1 

1043.7 

2.214 

19.12 

104 

2 

212.2 

1055.0 

2.257 

20.12 

105 

4 

215.0 

1067.0 

2.288 

21.12 

106 

2 

221.0 

1076.8 

2.352 

22.12 

107 

2 

226.5 

1086.9 

2.411 

23.12 

107 

8 

235.2 

1094.3 

2.503 

24.12 

108 

0 

249.8 

1097.1 

2.659 

25.12 

107 

5 

268.6 

1092.8 

2.858 

26.12 

106 

0 

294.2 

1078.4 

3.130 

27.12 

103 

9 

323.4 

1057.7 

3.438 

28,12 

101 

3 

355.7 

1032.7 

3.778 

29.12 

986.4 

388.2 

1005.9 

4.120 

30.12 

957.0 

423.6 

976.9 

4.490 

31.12 

926.2 

463.2 

946.4 

4.902 

32.12 

895.2 

506.5 

915.7 

5.351 

33.12 

863.6 

558.3 

884.4 

5.885 

34.12 

831, 

,6 

625.0 

852.6 

6.571 

35.12 

800.6 

709.7 

821.6 

7.439 

36.12 

772,3 

821.0 

793.3 

8.577 
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TABLE  III 


PREDICTION  OF  TRANSITIONAL  AND  TURBULENT 
BOUNDARY  LAYER  DEVELOPMENT  AND  HEAT  TRANSFER  RATES 


(Adjusted  Transition  Point,  S=8.4") 


n=2 

Experimental 
FLow  Field 


s 

NR0 

Nxt 

Nu 

9.4 

84J.6 

2.389 

10.4 

869.9 

2.245 

11.4 

890.5 

2.190 

12.4 

908.7 

2.162 

13.4 

926.8 

2.132 

14.4 

943.2 

2.119 

15.4 

959.2 

2.109 

16.4 

974.2 

2.109 

17.4 

988.4 

2.116 

18.4 

1000.3 

2,146 

19.4 

1011.6 

2.185 

20.4 

1022.6 

2.227 

21.4 

1032.2 

2.288 

22.4 

1041.7 

2.355 

23.4 

1047.7 

2.463 

24.4 

1049.2 

2.625 

25.4 

1042.5 

2.837 

26.4 

1026.9 

3.105 

27.4 

1005.9 

3.406 

28.4 

981.5 

3.725 

29.4 

955.4 

4.050 

30.4 

927.1 

4.402 

31.4 

897.7 

4.737 

32.4 

868.0 

5.214 

33.4 

837.4 

5.734 

34.4 

806.7 

6.391 

35.4 

776.4 

7.247 

36.4 

750.8 

8.234 
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TABLE  IV 


EXPERIMENTAL  MEASUREMENTS  AND  RESULTS 


For  the  Test  Point: 

M  =  8.0  P  =  590  psia 

00  boo 

T„  =  1 884°R,  Tw  =  540°R 

aoo 


S 

PxlO3 

2.615 

0.4471 

3.615 

0.4643 

4.115 

0.4541 

5.115 

0.3075 

5.615 

0.1754 

6.115 

0.0688 

7.365 

0.0344 

8.865 

0.0193 

10.865 

0.0206 

12.865 

0.0241 

14.865 

0.0392 

16.865 

0.0268 

18.865 

0.0227 

19.865 

0.0172 

24.365 

0.0235 

25.015 

0.0265 

26.615 

0.0330 

27.365 

0.0357 

28.115 

0.0380 

28.865 

0,0410 

30.365 

0,0460 

31.865 

0.0520 

33.365 

0.0610 

34.115 

0.0670 

34.865 

0.0760 

N^xl  0 

NNu 

146.6 

227.5 

209.9 

260.5 

243.7 

- 

205.5 

309.0 

136.0 

197.5 

62.5 

96.5 

33.6 

54.6 

41.8 

46.2 

51.0 

62.9 

72.0 

103.5 

67.6 

119.0 

78.2 

155.0 

88.0 

191.0 

94.0 

230.0 

102.4 

173.2 

117.8 

361.0 

145.3 

457.0 

189.0 

545.0 

208.2 

640.0 

221 J. 

742.0 

263.9 

894.0 

433.0 

1115.0 

480.4 

1334.0 

658.9 

1417.0 

74-1.2 

- 
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FIG.  lb.  The  Cone-Cylinder-Flare  Configuration,  Conical  Nose  Portion  and  Shoulde 
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FIG.  2.  Analytical  and  Experimental  Pressure  Distributions 


CONE  1 
CYLINDER 
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versus  S(inches).  Test  1. 


40 


FIG.  3b.  Heat  Transfer  Distribution  Over  Body  in  Terms  of  2  versus  S(inches). 

Test  2. 
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FIG.  4.  Heat  Transfer  Distribution  Over  Body  in  Terms  ar  2  versus  S  (inches).  Test  3 


S  INCHES 

FIG  5.  Heat  Transfer  Distribution  Over  Body  in  Terms  of  versus 

S (inches).  Test  4. 
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4.0 
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FIG.  7.  Heat  Transfer  Results  Over  Body  In  Terms  of  Nu(  versus  N* ,  Test  2 
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FIG.  8.  Heat  Transfer  Results  Over  Body  in  Terms  of  Nu'  versus  N'  ,  Test  3 
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FIG.  9.  Heat  Transfer  Results  Over  Body  in  Terms  of  Nu'  versus  N'  ,  Test  4 
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TRANSFORMED  MOMENTUM  THICKNESS,  ® 


SIMILARITY  PARAMETER,  £ 

FIG,  12.  The  Transformed  Momentum  Thickness  at  an 
Enthalpy  Ratio  gw=0.2772 
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SURFACE  DISTANCE  S 

FIG.  1 3a.  Distribution  of  the  Reynolds  Number  N  .in  the  Laminar  Region 
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SURFACE  DISTANCE,  S 

FIG.  13b.  The  OveraLL  Distribution  of  the  Reynolds  Number  N 
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FIG.  14.  Piecewise  Distribution  of  ~B  Over  the  Cone -Cylinder  Shoulder 


FREE  STREAM  STAGNATION  PRESSURE,  PSoo  PSIA. 
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TRANSITION  POINT,  §TRANs 

FIG*  15,  The  Effect  of  Stagnation  Pressure  on  the  Laminar  to 
Turbulent  Transition  Point 
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The  model  chosen  for  this  investigation  is  a  cone -cylinder- 
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portion  smoothly  blended  by  a  shoulder  radius  into  a  long  cylindrical 
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The  model  was  rested  at  a  free  stream  Mach  number  of  8  and  a 
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obtained  by  cooling  the  model  prior  to  the  test  with  liquid  nitrogen. 

The  stagnation -to -wall  temperature  ratios  ranged  from  10  to  3,3, 

The  theoretical  predictions  give  good  results  for  the  heat 
transfer  rates  in  the  laminar  region,  and  fair  prediction  in  the 
transitional  and  turbulent  regimes  extending  over  the  shoulder  and 
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correct,  and  underestimate  the  heating  rate  by  a  factor  as  high  as  2, 
Conversely,  the  "fiat  pia’e  reference  enthalpy"  method  is  found  more 
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increasingly  to  overestimate  the  heating  rates  over  the  forward  portion 
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